Design and Simulation of THz Quantum Cascade Lasers 
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Strategies and concepts for the design of THz emitters based on the quantum cascade scheme 
are analyzed and modeled in terms of a fully three-dimensional Monte Carlo approach; this allows 
for a proper inclusion of both carrier-carrier and carrier-phonon scattering mechanisms. Starting 
from the simulation of previously published far-infrared emitters, where no population inversion is 
achieved, two innovative designs are proposed. The first one follows the well-established chirped- 
superlattice scheme whereas the second one employs a double-quantum well superlattice to allow 
energy relaxation through optical phonon emission. For both cases a significant population inversion 
is predicted at temperatures up to 80 K. 



Since the first demonstration of the Quantum Cascade 
(QC) laser in 1994|l|, its performance has experienced 
tremendous improvement , and the range of emission 
wavelengths has been continuously extended^, ||, ^, 0- 
At present, the longest wavelength QC lasers operate at 
A ~ 21.5 /im and 24 /mi[j7|, still above the LO-phonon 
energy threshold of the host material (Tilolo ~ 34 meV 
in InGaAs, fiiOLO ~ 36 meV in GaAs). Although lasing 
at longer wavelength has not been observed yet, electrolu- 
minescence in the THz region of the spectrum has already 
been detected from a variety of QC structu res |[ ||, [Io| as 
well as from other quantum well device s %A\ |l2|, stimulat- 
ing a number of theoretical proposals [p. 3|. These strong 
efforts aimed at fabricating a THz semiconductor laser 
are mainly driven by the lack of compact, convenient 
solid-state sources operating at THz frequencies, despite 
the many possible applications in wireless communica- 
tions, medical imaging, security screening, etc. 

As prototypical design for THz emitters, we shall 
consider the GaAs/AlGaAs QC structure proposed by 
Rochat et al.^. This can be regarded as a scaled 
down version of the conventional mid-infrared QC de- 
sign, based on a vertical optical transition (see Fig. 1 
in Ref. || for its band diagram). This implementation is 
characterized by a narrow injector miniband, chosen to 
avoid scattering of electrons from the upper state through 
LO-phonon emission and cross-absorption of the emitted 
light. Moreover, it features a relatively small electron 
tunneling between the states in the active region and in 
the injector. As stressed before, while electroluminescent 
devices were reported by several groups H [l^], to date no 
population inversion could be achieved. In fact, several 
problems inherent to this design undermine the function- 
ality of the structure. First, the width of the miniband 



forbids the use of LO-phonon emission as the main mech- 
anism to deplete the lower state of the lasing transition. 
Second, the small miniband dispersion and low tunneling 
probability also restrict efficient electrical transport to a 
narrow range of currents and voltages, limited by the on- 
set of negative differential resistance. Third, the small 
energy separation between the levels in the active region 
and the quasi-Fermi level of the injector electrons makes 
the device very sensitive to hot-carrier effects ('backfill- 
ing' of the upstream active region). 

In order to gain a better understanding of these limit- 
ing factors, we have applied a g lobal simulation scheme, 
recently proposed in Refs. [lj, [D| to the prototypical 
structure of Rochat et aL||. Such a technique consists 
of a Monte Carlo sampling of a coupled set of fully three 
dimensional (3D) Boltzmann-like equations. To properly 
model energy relaxation as well as carrier thermaliza- 
tion, we included both carrier-phonon (c-p) and carrier- 
carrier (c-c) interaction mechanisms. Thanks to peri- 
odic boundary conditions, we are able to get the device 
current-voltage characteristics as well as its gain spec- 
trum without resorting to phenomenological parameters. 

At 80 K, under the design electric field of 2.7 kV/cm 
and assuming a sheet density of 1.2 xlO 10 cm~ 2 , our 
simulated experiments give an electron density in the 
upper laser level 2 of about 8.5 x 10 s cm~ 2 , compared 
to the 3.5 x 10 9 cm -2 in the lower state 1. These values 
are far from the transparency condition of the intersub- 
band transition, as already inferred from experimental 
data. This is not a consequence of rapid non-radiative 
depopulation of state 2 but rather of a slow extraction 
of electrons from level 1 into the downstream injector re- 
gion. Indeed our simulation shows rather a long lifetime 
p6[ of level 2, — 2.2 ps, mainly determined by c-p 
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scattering of electrons with high in-plane momentum to 
subband 1 (r 21 = 2.7 ps). On the other hand for level 

1 we compute a significantly longer lifetime t\ = 14.5 ps 
resulting from almost comparable c-c and c-p contribu- 
tions. In reality, the electron dynamics is dominated by 
very fast c-c scattering with a quasi-resonant level (1' in 
Ref. 7). However, such a process acts in both directions 
(1 — > 1' and 1' — > 1) with nearly equal rates, thus giving 
only a marginal contribution to the electron depletion of 
state 1. From the simulation we also obtain the value of 
the operational current density, 50A/cm 2 , in very good 
agreement with the measured one at this bias field. 

Starting from these considerations we have designed 
two structures tuned for emission at A ~ 69 fim that can 
overcome the limitations associated with the above pro- 
totypical design. Both structures are based on a vertical- 
transition configuration, which is known to lead to larger 
dipole matrix elements and narrower linewidths. The 
first scheme employs a conventional chirped-superlatticc 
design|l7]|, which allows flat minibands in the active re- 
gion without requiring large concentrations of dopants. 
The band diagram under an average applied electric field 
of 3.5 kV/cm is shown in Fig. ffl. Its operating strat- 
egy is based on the same concepts successfully imple- 
mented for A ~ 17 — 24 /jm QC lasers j| ^,0. In order 
to minimize the density of electrons in the lower laser 
subband 1 we employ a dense miniband with seven sub- 
bands, which provides a large phase space where electrons 
scattered either from subband 2 or directly from the in- 
jector can spread. The miniband dispersion is chosen 
as large as possible compatibly with the need of avoid- 
ing cross-absorption. This suppresses thermal backfilling, 
and provides a large operating range of currents and volt- 
ages. Also in this structure, however, energy relaxation 
within the first miniband appears to be hindered by the 
lack of final states with appropriate energy to allow for 
LO-phonon emission. Nevertheless, carrier-carrier inter- 
actions may beneficially operate as an activation mech- 
anism, which can provide sufficient in-plane momentum 
for the electrons to open a scattering path via optical 
phonon to the lower subbands. This can be clearly seen 
from Tab. 1, where we report the simulated steady-state 
distribution of carriers at 80 K in the various subbands, 
with and without c-c scattering in the simulation. As ex- 
pected, in the absence of c-c interaction, the carrier den- 
sity remains high in the upper states of the injector mini- 
band and particularly in subband 1. On the other hand, 
in the presence of c-c scattering, electrons efficiently re- 
lax into the lower states of the injector (g,A), from where 
they are readily transferred into the upper laser state of 
the following period. In the latter case a population in- 
version of An = 1.1 X 10 9 cm~ 2 results between subbands 

2 and 1: compared to the previous case, the number of 
electrons of subband 1 is decreased by 40% whereas that 
of state 2 increases by a factor of four. Calculated life- 
times are r 2 = 0.8 ps, t 2 i = 8.3 ps, t\ — 2.2 psjl^] with r 2 
and r 2 i dominated by c-p scattering processes. Note that 
t 2 i is much greater than t 2 due to the superlattice nature 



of the active region. This implies that t\ < t 2 i, which 
is the necessary condition for population inversion. Fur- 
thermore, the number of carriers escaping via subband 3 
is negligible. 

In order to further improve the extraction rate from 
the lower laser state 1, we designed another structure 
which allows for a direct use of c-p scattering as the de- 
pleting mechanism without compromising the lifetime of 
the upper state 2. This is achieved by using a double- 
quantum well superlattice (DQW-SL)|^9|, which leads to 
a split miniband in the injector. By a careful design of 
the DQW-SL it is possible to separately adjust the sub- 
minibands dispersion and separation, where the latter 
can be chosen to match the LO phonon energy. Thereby 
we avoid cross- absorption and, at the same time, recover 
the possibility of using the efficient direct scattering by 
LO-phonons as the main relaxation process. To this end, 
we have also tried to delocalize the lower laser state 1 
over both the active region and the following injector. 
Although this leads to a slight reduction of the dipole 
matrix element with the upper laser state, it increases 
significantly the LO-phonon assisted transitions to the 
injector subbands, in particular A, B, g. In that way 
the lifetime of state 1 is kept as short as 0.3 ps, resulting 
from both c-c and c-p scattering. In order to preserve the 
long lifetime of the upper laser level 2, the latter is well 
confined in only two wells and is spatially separated from 
the injector. This keeps the scattering probability of elec- 
trons out of level 2 rather low, leading to a relatively long 
lifetime of 1.1 ps (r 2 i = 1.9 ps), which is dominated by 
c-p scattering. As reported in Tab. 2, population inver- 
sion in this structure is achieved even without taking into 
account c-c scattering. This confirms the good function- 
ality of the phonon-split sub-minibands. In the presence 
of c-c scattering, the population of the lower laser level 
increases due to faster scattering from level 2 to 1. How- 
ever, this is compensated by a stronger injection from 
subband g into subband 2, resulting in an increased pop- 
ulation of the latter. The population inversion amounts 
to about 10 percent of the total sheet density per period. 
Note that, owing to the larger voltage drop per period, 
thermal 'backfilling' of carriers is only a minor problem, 
thereby allowing for a higher doping in this second de- 
sign. 

Considering the dipole matrix elements of 7.8 nm 
and 4.5 nm computed for the optical transitions 
in the two structures, respectively, and assuming a 
linewidth of 2 meV as experimentally detected at sim- 
ilar wavelengths ||, we estimate, on the basis of the 
above populations, gain coefficients of about 31 cm -1 and 
110 cm -1 . These values compare favorably to the opti- 
cal losses (approximately 50 cm -1 ) measured in double- 
surface plasmon waveguides at THz frequencies |2C| , and 
indicate that laser action with such active region designs 
is a realistic goal. 

In summary, we have proposed two QC structures 
designed for THz lasing. Our theoretical analysis of 
their performance is based on a global simulation scheme 
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which allows for a proper inclusion of the most relevant 
scattering mechanisms. The proposed designs both ex- 
hibit significant population inversion up to 80 K. This, 
together with the large transition dipolcs, is expected to 
lead to high optical gain. Finally, it should be noted that 
our simulation predicts current densities of more than 1 
kA/cm 2 which is roughly one order of magnitude higher 
than the current density measured in the prototypical 
design of Ref. before the onset of negative differen- 



tial resistance. We attribute this effect to the improved 
carrier relaxation and higher tunneling efficiency of our 
structures. 

This work was supported in part by the European 
Commission through the FET project WANTED. R. K. 
acknowledges support by the C.N.R. The authors would 
like to thank Fabio Beltram for helpful discussions and a 
careful reading of the manuscript. 



[1] J. Faist, F. Capasso, D. L. Sivco, C. Sirtori, A. L. 

Hutchinson, A. Y. Cho, Science 264, 553 (1994). 
[2] F. Capasso, C. Gmachl, R. Paiella, A. Tredicucci, A. L. 

Hutchinson, D. L. Sivco, J. N. Baillargeon, A. Y. Cho, 

H. C. Liu, IEEE J. Select. Topics Quantum Electron. 6, 

931 (2000). 

[3] D. Hofstaetter, M. Beck, T. Aellen, J. Faist, U. Oesterle, 
M. Illegems, E. Gini, H. Melchior, Appl. Phys. Lett. 78, 
1964 (2001). 

[4] J. Faist, F, Capasso, D. L. Sivco, S. N. G. Chu, A. Y. 

Cho, Appl. Phys. Lett. 72, 680 (1998). 
[5] A. Tredicucci, C. Gmachl, F. Capasso, D. L. Sivco, A. L. 

Hutchinson, A. Y. Cho, Appl. Phys. Lett. 74, 638 (1999). 
[6] A. Fredicucci, C. Gmachl, M. C. Wanke, F. Capasso, A. 
L. Hutchinson, D. L. Sivco, S. N. G. Chu, A. Y. Cho, 
Appl. Phys. Lett. 77, 2286 (2000). 
[7] R. Colombelli, F. Capasso, C. Gmachl, A. L. Hutchinson, 
D. L. Sivco, A. Tredicucci, M. C. Wanke, A. M. Sergent, 
A. Y. Cho, Appl. Phys. Lett. 78, 2620 (2001). 
[8] M. Rochat, J. Faist, M. Beck, U. Oesterle, M. Ilegems, 

Appl. Phys. Lett. 73, 3724 (1998). 
[9] B. S. Williams, B. Xu, Q. Hu, M. R. Melloch, Appl. Phys. 
Lett. 75, 2927 (1999). 

[10] J. Ulrich, R. Zobl, K. Unterrainer, G. Strasser, E. Gornik, 
Appl. Phys. Lett. 76, 19 (2000). 

[11] W. Hilber, M Helm, K Alavi, R. N. Pathak, Appl. Phys. 
Lett. 69, 2528 (1996). 

[12] P. Bakshi, K. Kempa, A. Scorupsky, C. G. Du, G. Feng, 
R. Zobl, G. Strasser, C. Rauch, Ch. Pacher, K. Unter- 
rainer, E. Gornik, Appl. Phys. Lett. 75, 1685 (1999). 

[13] P. Harrison, R. A. Soref, IEEE J. Quantum Electron. 
QE-37, 152 (2001) and references therein. 

[14] R. C. Iotti, F. Rossi, Appl. Phys. Lett. 78, 2902 (2001) 

[15] R. C. Iotti, F. Rossi, Phys. Rev. Lett., in press. 

[16] Effective inter-level scattering rates and lifetimes can be 
directly sampled from our fully 3D simulated experi- 
ments. In this way we can also properly account for the 
strong balance between in- and out- scattering typical 
of quasi-elastic interaction processes. This allows us to 
translate our microscopic description into the more con- 
ventional n-level rate-equation language. 

[17] A. Tredicucci, F. Capasso, C. Gmachl, D. L. Sivco, A. 
L. Hutchinson, A. Y. Cho, Appl. Phys. Lett. 73, 2101 
(1998) 

[18] As for the structure of Rochat et al. also in this case there 
is fast exchange of electrons with a third quasi-resonant 



state. 

[19] M. C. Wanke, F. Capasso, C. Gmachl, A. Tredicucci, D. 
L. Sivco, A. L. Hutchinson, S. N. G. Chu, A. Y. Cho, 
IEEE Photon. Technol. Lett. 13, 278 (2001). 

TABLE I: Population of the individual levels in the 'chirped' 
SL, calculated at 80 K with only carrier-phonon (c-p) scat- 
tering and with both c-p and carrier-carrier (c-c) scattering. 
Letters refer to the injector states comprising the first mini- 
band in Fig. |l|. g represents the injector ground state while 
capital letters label the other levels in ascending order. 



Subband 
index 


Energy 
(meV) 


Pop. (c-p only) 
(10 9 cm- 2 ) 


Pop. (c-p & c-c) 
(10 9 cm" 2 ) 


3 


51.7 


0.02 


0.06 


2 


36.1 


0.94 


3.99 


1 


17.9 


4.72 


2.83 


F 


14.1 


7.13 


4.72 


E 


11.0 


8.08 


6.08 


D 


7.7 


7.34 


5.45 


C 


4.5 


5.14 


5.98 


B 


3.2 


1.99 


4.62 


g 





6.19 


7.87 



TABLE II: Population of the individual levels in the DQW 
SL, calculated at 80 K with only carrier-phonon (c-p) scat- 
tering and with both c-p and carrier-carrier (c-c) scattering. 
Subbands are named with the same convention as in table fl 



Subband 
index 


Energy 
(meV) 


Pop. (c-p only) 
(10 10 cm- 2 ) 


Pop. (c-p & c-c) 
(10 10 cm- 2 ) 


2 


63.4 


0.27 


1.22 


1 


45.5 


0.07 


0.46 


E 


42.6 


0.02 


0.27 


D 


37.1 


0.05 


0.27 


C 


28.7 


0.27 


0.27 


B 


9.2 


1.22 


1.16 


A 


5.9 


0.82 


1.43 


g 





5.24 


2.99 
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FIG. 1: Conduction band energy diagram of the 'chirped' 
superlattice under an average electric field of 3.5 kV/cm. The 
layer thickness (in nm) are, from left to right, starting from 
the injection barrier 4.3/ 18.8/ 0.8/ 15.8/ 0.6/ 11.7/ 2.5/ 
10.3/ 2.9/ 10.2/ 3.0/ 10.8/ 3.3/ 9.9, where Alo.15Gao.85 As 
layers are in bold face and the 10.2 nm wide GaAs well is 
doped 4 x 10 16 cm~ 3 . Also shown are the moduli squared 
of the wavefunctions; the optical transition occurs between 
states 2 and 1. 




FIG. 2: Conduction band energy diagram of the DQW- 
SL under an average electric field of 8.7 kV/cm. The layer 
thickness (in nm) are, from left to right, starting from the 
injection barrier 3.5/ 6.8/ 1.7/ 6.4/ 2.5/ 7.9/ 0.6/ 7.5/ 2.0/ 
6.8/ 1.0/ 6J>/ 2.0/ 5.9, 1.4/ 5.5, where Alo.3Gao.7As layers 
are in bold face and doped (6 x 10 16 cm -3 ) GaAs layers are 
underlined. The Al concentration in this design, and thus the 
barrier height, is raised compared to the previous design in 
order to better accomodate the split miniband. 



